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Abstract. Quantum wire networks have recently become of great interest. 
Here we deal with a novel nano material structure of a Double Gyroid wire 
network. We use methods of commutative and non-commutative geometry 
to describe this wire network. Its non-commutative geometry is closely 
related to non-commutative 3-tori as we discuss in detail. 



Introduction 

Interfaces that can be modeled by surfaces of constant mean curvature (CMC) 
are ubiquitous in nature and can now be synthesized in laboratory. Recently, 
Urade et al. [23] have reported fabrication of a nano-porous silica film whose 
structure is related to a specific CMC surface — the Gyroid. The structure is 
three-dimensionally periodic and has three components}^ a thick surface and 
two channels, as detailed below. The interface between the wall and the chan- 
nels approximates a Double Gyroid. 

Urade et al. [23] have also demonstrated a nanofabrication technique in which 
the channels are filled with a metal, while the silica wall can be either left in 
place or removed. These novel materials open a wide field of applications due 
to their topological and geometric features. The channels are a few nanome- 
ters wide and, when filled with a conducting or semiconducting material, are 
expected to acquire certain characteristics of one-dimensional quantum wires 
(such as a blueshift of the spectrum and an enhanced density of states), while 
remaining three-dimensional in other respects |13j . Geometrically, the one- 
dimensional structure appears since each of the channels can be retracted to a 
skeletal graph [lOj, the Gyroid graph. 

We will concentrate on the resulting geometry and topology of these networks 
in the possible presence of a constant magnetic field. Topological quantities are 
of particular physical interest as they remain stable under continuous pertur- 
bations. In practice, those should be small as not to break the structure. 

We use two approaches: the first is purely classical and the second is the 
quantum/non-commutative approach due to Alain Connes [S]. The classical 
results we provide are a study of the fundamental group of the channel. Here 
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we determine the fundamental group of each channel system. This group is 
the commutator subgroup of the free group in three generators. A surprising 
result coming directly from the Gyroid geometry is that there is a new length 
function on the free group which is different from the ordinary word length. 
With the help of this length function we determine the shortest loops at any 
given point. There are 15 such topologically distinct shortest loops (30 if one 
includes orientation). They split into two groups which are distinguished by 
their cyclic symmetry which is either of order two or three. In both groups 
there are three generators using this additional symmetry. These loops are of 
particular interest since numerical simulations |13] show the possibility of an 
enhanced density of states in a double-gyroid quantum wire, due to states that 
are nearly localized near such loops. We also calculate the flux of a constant 
magnetic field through these loops. The tool is an effective unit vector. The 
result of the calculation is that these effective unit vectors have a particularly 
simple form (see Table [T]) . This fact should also be relevant for the study of the 
spin-orbit coupling of the loop-localized states. 

Our study of the non-commutative geometry is motivated by one of the big 
early successes of the non-commutative geometry of Alain Connes. This was 
the description by Bellissard et al. of the quantum Hall effect [3J. This al- 
lowed to explain the integer effect in terms of the non-commutative geometric 
properties. The underlying geometry in that situation is the quantum 2-torus. 
Recently there have been further analyses on the fractional effect using hyper- 
bolic geometry as a model [15] . The conceptual approach as outlined in [2] is to 
replace the Brillouin zone by a non-commutative Brillouin zone which is given 
by a C* -algebra that contains the translational symmetry operators and the 
Hamiltonian. In this geometry relevant quantities can often be expressed in 
terms of the K theory of this algebra. This Abelian group captures informa- 
tion related to the the topology or better the homotopy type of the algebra or 
geometric setup. These are again quantities which are stable under continuous 
deformations. A prime example is the Hall conductance. 

A general fact which is pertinent to our discussion is that the ii'-groups also 
serve to label the gaps in the spectrum. Roughly this goes as follows. Given a 
gap in the spectrum there is a projector projecting to the energy levels below 
the gap. This projector in turn gives rise to a /C-theory element. If one knows 
the ordered i^-theory then one can also sometimes deduce if only finitely or 
infinitely many gaps are possible. 

In our situation, we determine the said C* algebra for one channel in the 
presence of a magnetic field and describe the iT-theories. The Hamiltonian 
we use is the generalization [20' '2] of the Harper Hamiltonian [12\ [2] adapted 
to our situatioij^ We call the resulting algebra the Bellissard-Harper algebra 
and denote it by This geometry is closely related to the non-commutative 
3-torus Tq. Here is a skew-symmetric 3x3 matrix determined by the 



The generalized Harper operator is also the operator underlying the non-commutative 
geometry of the quantum Hall effect [31 



THE GEOMETRY OF THE DOUBLE GYROID WIRE NETWORK 



3 



magnetic field. In fact we show that the algebra is isomorphic to a subalgebra 
of the 4x4 matrix algebra with coefficients in the non-commutative 3-torus. 
By varying the magnetic field, we obtain a three-parameter family of algebras. 
We prove that at all but finitely many points this algebra is the full matrix 
algebra and hence Morita equivalent to Tq. Since Tq is simple at irrational Q 
one would expect this to be true on a dense set (see §3). We not only prove 
that this expectation holds at the irrational points, but are able to extend 
this result to almost all rational points. At certain special values of the field 
which we enumerate, however, the algebra is genuinely smaller leading to a 
possibly different ET-theory. At these points the material may also exhibit 
special properties. 

The ordered K-theory of the 3-torus is completely known [181 [13 E] and 
actually completely classifies the isomorphism classes of such tori [19\ [7] . Our 
algebra always injects into a Tq. From this and the knowledge of the ordered 
iC-group of the non-commutative torus, we obtain the result that there are 
only finitely many gaps possible at rational values of 0. 

Our approach to both the classical and the quantum geometry uses graph 
theoretical methods. The relevant graphs are the Gyroid graph which is a 3- 
regular graph and its quotient under the translation symmetry group. In order 
to give the matrix algebras explicitly we introduce the notions of a graph rep- 
resentation. Using rooted spanning trees, we are able to represent the algebra 
^ in terms of matrices, where they are amenable to direct computations. 

As the geometry of the Gyroid and its channels is quite difficult, we also 
treat the honeycomb lattice as a two-dimensional analogue. Indeed, the Gyroid 
lattice graph (the graph onto which each channel retracts) is in many ways 
the three dimensional analogue of the honeycomb lattice, which is why when 
developing the more general parts of the theory we will parallely consider these 
two cases as our main examples. Both graphs are 3-regular, and both of the 
graphs are not mathematical lattices but only physical lattice^ This means 
that they give rise to two groups, one which is the space or symmetry group 
and the other is the group of the lattice of which they can be considered a 
subset. In the honeycomb case these are both but embedded into as the 
triangular lattice and its dual. For the Gyroid the groups are both Z^, but they 
are embedded as a body centered cubic (bcc) and a face centered cubic (fee) 
lattice in M.^ which are again dual to each other. Furthermore the fundamental 
group of the honeycomb lattice is the commutator subgroup of the free group 
on two elements, while for the Gyroid tti is the commutator subgroup of the 
free group in three generators. And in both cases, we find a new length function 
on the free group induced by the geometry of the lattice. 

The paper is organized as follows: 

In the first section, we start with a review of the classical geometry of the 
Gyroid and then prove the results on the fundamental group and the small- 
est loops. In the second section, we formally introduce graphs and lattices 



^See isii 



for the disambiguation of the use of the word "lattice" . 



4 



RALPH M. KAUFMANN, SERGEI KHLEBNIKOV, AND BIRGIT KAUFMANN 



and the relevant groups associated to them. We then define their representa- 
tion in Hilbert spaces and the associated Harper operators. Finahy, we show 
how to obtain a matrix representation of using rooted spanning trees. In 
the third section, we apply the general theory of the second section to the 
non-commutative geometry of a lattice using a Harper operator, detailing the 
honeycomb and the Gyroid case. Here we also briefly review projective repre- 
sentations and explain how they arise in the presence of a magnetic field. In 
this context, we can already show that there can only be finitely many gaps for 
the Gyroid at rational 0. We also outline the general approach to calculating 
^ and its JC-theory and the expected results. In particular, we compute the 
-fC-theory in the commutative case in terms of a cover of a torus. 

In the last section we apply the outlined strategies to calculate the algebras 
and their X-theory for Bravais lattices, the honeycomb lattice and the Gyroid 
lattice graph. 

1. The CLASSICAL GEOMETRY OF THE DOUBLE GYROID 

1.1. The double Gyroid and its channels. The Gyroid is an embedded 
CMC surface in [9j. It was discovered by Alan Schoen [21]. In nature it 
was observed as an interface for di-block co-polymers |ill] [22] . The interface 
actually consists of two disconnected surfaces. Each of them is a Gyroid surface. 
The Double Gyroid (DC) is a particular configuration of two mutually non- 
intersecting embedded Gyroids. 

A single Gyroid has symmetry group I4i32 while the double Gyroid has the 
symmetry group laSd where the extra symmetry comes from interchanging the 
two Gyroids]^ 

Since CMC surfaces are mathematically hard to handle level surfaces have 
been suggested as a possible approximation in [2]. The level surface approx- 
imation for the Double Gyroid again consists of two level surfaces. We will 
call the two surface interfaces S\ and S^- For the discussion at hand it is not 
relevant if the two surfaces are actually the CMC surfaces or their level surface 
approximations. One example of a DC approximation is given by the family of 
level surfaces [H] 

Lt : sin X cosy + sin y cos z + sin z cos x = t (1-1) 
A model for the double Gyroid is then given by and for < it; < \/2. 
The complement \ G of a single Gyroid G has two components. These 
components will be called the Gyroid wire systems or channels. 

There are two distinct channels, one left and one right handed. Each of these 
channels contracts onto a graph, called skeletal graph in j2H I10| . We will call 
these graphs r+ and r_. Each graph is periodic and trivalent. We fix F^ to 
be the graph which has the node vq = (f, fi §)• We will give more details on 
the graph F+ below. 



Here /4i32 and Ia3d are given in the international or Hermann-Mauguin notation for 
symmetry groups, see e.g. [8]. 
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lil 

Figure 1 . The Channel C+ 

The channel containing r+ is shown in Figure [TJ A (crystal) unit cell of the 
channel together with the embedded graph r+ is shown in Figure [2] and just 
the graph is contained in Figure [3j 

In constructing the surface, one can actually start with one of these graphs 
and evolve the surface from it [lOj. The symmetry of one skeletal graph is that 
of a single Gyroid /4i32. 

To obtain the double Gyroid one can evolve from both the graphs. The 
common symmetry group of both of these graphs is laSd. This symmetry 
group can actually be defined already on the nodes V± of the two graphs F-i-. 
The subgroup I4i32 is then determined to be the subgroup that fixes both sets 
of nodes (setwise). 

In the case of the double Gyroid S = SiJI S2, the complement C = 'R.^\S of 
the disjoint union of the two Gyroid surfaces has three connected components. 
These are two channel systems C+ and C_, each of which can be retracted to 
its skeletal graph F-i-. In fact, the skeletal graph is a deformation retract [lOj . 
There is a third connected component F and we will call F = FU S. This is a 
3-manifold with two boundary components, more precisely, dF = S = SiU 82- 
F can be thought of as a "thickened" (fat) surface. The thickness is fixed by 
the parameter w. Note that F can be retracted to each of the two boundary 
surfaces Si. In fact there is a deformation retract of F onto the Gyroid. 
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Figure 2. The graph r+ embedded into its channel C-^ 



This means that all homotopical information about the complements of the 
double Gyroid is encoded in the Gyroid/level surface and the two skeletal 
graphs. 

The translational symmetry group for both the Gyroid and the double Gyroid 
is the Bravais lattice bcc. Note that we will usually use "lattice" in the physical 
terminology, e.g. speak about the honeycomb lattice. The term "Bravais lattice" 
will be used to denote a maximal rank mathematical lattice, i.e. a free rank n 
Abelian subgroup of M". In order to preempt any confusion, we give precise 



definitions for our terminology in ^3.2 



1.2. The skeletal graph. 



1.2.1. The vertices and edges. We will now describe the graph r+ em- 
bedded into M^. 




Figure 3. The graph r+ 

Set 

^0 = (g,g,g) V, = {-,-,-) 

,3 7 5, ,17 3, 

-1 = (8,8,3) -5 = (8,8'8) 

3 17 11 1 

^2 = (3, 8,3) ^6 = (8,8'8) 

,537, -731n /-ir\\ 

-3 = (8,8'8) ^^ = ^8'8'8) ^'-'^ 
and let V+ = {vq, . . . , vy}. 

Furthermore acts on by translations and we let = Z^(y+) be the 
image of the set V+ under this action. We will sometimes call this set of points 
the Gyroid lattice. 

Given two points u, w G let vw = {(1 — t)v + tw\t G [0, 1]} be the line 
segment joining them. Also given a point v we let Tx{v) =v + {l, 0, 0),Ty{v) = 
V + (0, 1, 0),T^{v) = v + {0, 0, 1) be the translated points. 

Consider the following set of line segments 
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f 4 {V5), V7 Tx (vq) , Vl Ty {V2), 

V5Ty{ve),V2T,ive),V3T,{v7)} (1.3) 

Let E+ = Z^{E), where again c acts as a subgroup of the translation 
group. 

Definition 1.1. The graph r+ is the graph whose vertices are V^, whose edges 
are E^ with the obvious incidence relations. 

We recall: 

Proposition 1.2. [lOj C+ can be deformation retracted onto the graph and 
F-i- is the component of the critical graph contained in C+. 

Corollary 1.3. The homotopy type of the complement T = W^\G is the same 
as that of two copies o/r+. In particular each channel has the same homotopy 
type as r+. 

This implies that all topological invariants C± which are homotopy invariant 
are isomorphic to those of r+. In particular, this means all homology, homotopy 
and K-groups of the topological space C+ or C_ are determined on F-i-. 

1.3. The quotient graphs. Quotienting out by either the standard transla- 
tion group or the bcc lattice, we obtain the following two quotient graphs. 

1.3.1. Crystallographic Quotient Graph. Let f^fT^'^*"' be the graph r+/Z^ 
thought of as an abstract graph. This graph has a natural map embedding into 
the 3-torus R^/Z^. 

Proposition 1.4. Y'^y^'-"-'' ig a cube embedded into the 3-torus. More precisely 
the vertices of T^^^'^*'"'^ are V and the edges are the images of E in M'^/Z^. 

The abstract graph is given in Figure |4] which also contains the images of the 
vectors Cj which name and orient all edges. 

Proof. Since all the vertices in the fundamental domain do not lie on the bound- 
ary, they give exactly the representatives of V in M^/Z'^. The classes of the edges 
Ey, are given exactly by the images of the set E in M^/Z^. □ 

1.3.2. The (maximal) Quotient Graph. By modding out by Z^, we have 
not yet used the full translational symmetry of F_|_ which is the bcc lattice. A 
set of generators for the bcc lattice is 

/i:= (1,0,0), /2 = (0,1,0), /3:= ^(1,1,1) (1.4) 
another set of generators which is more symmetric and we will use later on is: 
51 = J(l,-l,l), 52 = J(-l,l,l), 53 = ^(1,1,-1) (1.5) 
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Figure 4. The graph r!|_''^**''' for the Gyroid 



C 






Figure 5. The complete square, the rooted spanning tree r 
(root A and edges 61,62,63) and the vectors corresponding to 
the oriented edges, the collapsed tree f+/r 



We let L = L{r+) be the free Abelian subgroup of generated by these 
vectors. We define r+ to be the abstract quotient graph T^/L. 

Proposition 1.5. is the graph with 4 vertices and 6 edges, where all pairs 
of distinct vertices are connected by exactly one edge. 

This graph is sometimes also called the complete square. Its incidence matrix 
has entries one everywhere except on the diagonal, where the entries are zero. 
This graph is shown in Figure [5j 

Proof. We see that there is an embedding of Z'^ C L, where L is the bcc lattice, 
so that we only have to divide f^'^y^^"'^ by the additional symmetry generated 
by d := Now mod Z^, = id and (the translation by d) simply 
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interchanges the vertices of Y'^^^''"''' as follows: o ug, ui -H- uy, U2 -H- V4 and 
V3 -H- t's- Hence we are left with 4 vertices, and we can choose the representatives 



vo, ■ ■ ■ ,V3. Checking the list of edges ( 1.3 ) one sees that indeed the 12 edges form 



pairs and one can choose the representatives ViVj; i ^ j : i, j = 0, . . . ,3. □ 

Notice that this corresponds to a Z/2Z symmetry of the graph p^'"?'''*"'^ n 
given by mapping each vertex to its diagonally opposite vertex and maps the 
edges accordingly. 

1.4. The Underlying Group and Lattice. There is another Bravais lattice 
hidden in the geometry of the Gyroid. This is the fee lattice. The nearest 
neighbor positions differ by vectors generating an fee lattice. 

This means in particular that after shifting by vq the positions of the vertices 
of r+ all lie on an fee lattice. 

In order to fix notation set: 

^0 e2 = V3- Vo 63 = ■U4 - Vo (1.6) 

V2 e5 =V7 - VA 66 = 1^5 - Vl (1.7) 

be the group of that is generated by 

62 = i I -iV "3 = J 1^ 0^ 

es = I I -1 I , ee = \(^0^ (1.9) 

(1.10) 

Proposition 1.7. The group r(r_|_) is isomorphic to 1? . The Bravais lattice 
it generates is a face centered cubic (fee). 

Proof. There are relations among the given as follows. 

ei = -e5 + e6, 62 = -64 - ee, 63 = 64 + 65 (l-H) 

so that we see that 63,64,65 generate. These vectors are linearly independent 
over Z since they are independent over M. Hence they provide a free basis 
and an isomorphism to Z^. The vectors 64, —65, —66 are the standard primitive 
vectors for the face centered cubic. □ 

Of course there are many other choices of basis here, e.g. {62, 64, 66}. 

Proposition 1.8. The vertices o/F-i- translated by —vq lie on the 3-dimensional 
face centered cubic lattice generated by 64,65,66. 

Proof. r+ is path connected and we take v^ as the base point. Each line seg- 
ments in E corresponds to an edge. Choosing an orientation for this edge 
defines a vector. Now, the statement follows from the fact that the vectors 
corresponding to the line segments in E and hence those of E'r^ are exactly the 
vectors ±61, ±66. □ 
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1.5. Fundamental Groups, Loops and Effective Normal Vectors. There 
is certain geometric information which can even already be read off from the 
simple graph r+. Such as its fundamental group or the minimal loops starting 
at a given vertex. A loop on r_|_ is minimal if it passes through a minimal 
number of edges. A more general treatment will be given in ^ 

Proposition 1.9. LetY^ be the free group in three variables. The (realization) 
of the graphs T-i- and r_|_ have following fundamental groups: 

(1) 7ri(r+) = [F3,F3] 

(2) 7ri(r+)=F3 

in particular r+ is the maximal Abelian cover o/r+. 

Since r_|_ is homotopic to one channel these results hold for each channel of 
the Gyroid. 

Proof. We start with r^. This graph is homotopic to the wedge product of 
three S^'s, whence the second claim follows. 



In view of Lemma 1.10 below (more generally by Proposition 3.5 ), we see that 
7ri(r_|_) is the subgroup that is the kernel of the map F3 to its Abelianization. 
Indeed the sum of powers of each of the generators in a word in the group 
has to be zero, which precisely means that such a word is in the kernel of the 
Abelianization map or in other words, the commutator group. □ 

Lemma 1.10. A loop on the graph r+ lifts to a loop on r_|_ if and only if each 
edge is traversed the same number of times in each direction. 

Proof. The "if" direction is clear since this means that the translations have to 



add up to zero. The equivalence follows from the general Proposition 3.5, by 
noticing that indeed the lift s h = 626^^6^^, I2 = 6165^63^,^3 = 636463 give 



rise to the vectors /j = /j of (1.5) and are linearly independent. □ 



Proposition 1.11. There are closed loops in the graph r_|_. Each minimal loop 
goes through 10 sites and at each point there are 30 oriented minimal loops or 
15 such undirected loops. 

Proof. A path in which one goes back and forth through an edge is homotopic 
to the path in which this step is omitted. By direct calculation one can see 
that such a path is given by traversing any five of the six edges. At the first 
step one has 3 choices of edges, at the second step there are two and then again 
two possibilities. Here one either returns to the original vertex or not. In the 
first case there are 2 completions and in the second case 3 completions to a 
minimal loop. Thus we have 2 • 3 • 5 = 30 possibilities. A detailed version of 



this enumeration is provided at the end of this section in 5.1, □ 



1.5.1. Explicit Calculation of the loops. Here we give the details of the 



calculation of the loops. Using Proposition 3.5 we have to look for paths that 
traverse each edge the same number of times. We first look at the cases where 
each edge is traversed only once in each direction. We call such a path a good 
path. We also have to keep the number of these edges minimal. This already 
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puts a simple constraint on the path. We may not go back and forth through 
the same edge. Starting at the given vertex we have to pass through 2 distinct 
edges. After this we have a choice, wc can cither go back to the starting vertex 
(case I) or we can go to the only vertex which we have not reached yet (case 
II). In case I the next oriented edge is fixed, but then we have two choices la 
and lb. After this wc have already used 5 edges so that the minimal possible 
number of oriented edges and hence the length of the loop is 10. Indeed, we 
can complete the edge path uniquely to a good path without traversing the 6th 
edge. 

In case II the fifth oriented edge for a good path is fixed, going opposite the 
first oriented edge. There is a choice for the sixth oriented edge in the path. 
Either returning to Vq (Hb) or not (Ila). The case lib has a unique choice for 
a fifths oriented edge for a good path and this has a unique completion to a 
good path involving 5 edges, again giving a length of 10 loops. In case Ila, we 
again have two choices for the fifth oriented edge (Hal) and (IIa2). Both these 
choices have a unique completions to good paths again of length 10. 

It remains to treat the cases where an edge is traversed more that once in 
each direction. Since we are not allowed to go back and forth on one edge the 
choices for the first three oriented edges are as above. Now in Case I we could 
go along the first oriented edge again, but this would lead us to traverse at least 
6 edges and hence would not be minimal. At the next step of case I the choices 
are precisely la or lb and traversing an edge twice in the same direction would 
not be minimal. For case II, the first stage where one could use an oriented 
edge twice for Ila is the fifth edge. But then one would need at least 6 edges 
counting multiplicities. After fixing the fifth oriented edge one can only increase 
the number of traversed oriented edges by not choosing a good path. Finally 
in case lib again the first edge with a choice to traverse an edge twice in one 
orientation is at the fifth oriented edge, but as before this would lead to a path 
of length greater than 10. 

So all in all we have 2-3 = 6 choices for the first two edges and then once 
these are fixed 5 choices for a good path. This means in all there are 30 such 
oriented paths. 

1.5.2. Explicit loops. On the graph r+ there is a symmetry group of order 
6 preserving the base point and the spanning tree r which permutes the vertices 
B,C and D. This is precisely the group that gives us the 6 first choices. 

In order to write down a shorter list, we will make the following observations. 
Since the inverse of a minimal path is a minimal path, we can cut down the 
number to 15. Now since each of the paths traverses five edges, it misses one. 
There are two cases: (1) the edge that is missed is incident to wq, i.e. 61,62 or 
63 and case (2) it is not, i.e. 63,64,60. Case (1) corresponds to Hal and IIa2, 
which case (2) corresponds to la, lb and lib. In the case (1) the vertex vq is 
traversed 1 additional time except at the start and end of the path and in case 
(2) it is traversed an additional 2 times. This decomposes the loops into two 
pieces of length 5 in case (1) or three pieces of length 3,3,4; 3,4,3; 4,3,3 in case 
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lib 



Figure 6. The combinatorial cases for the minimal loops of the Gyroid 

(2) for la, lb, lib respectively. With each loop, its cyclic permutation of these 
components is also a loop. There are 2 such loops in case (1) and 3 such loops 
in case (2). This also explains the 15 loops as 15 = 3 • 2 + 3 • 3 and permutes 
the cases Hal, IIa2; and Ia,IIb and lb cyclically. 

We give here a list of six basic loops which may be useful for further dis- 
cussion. To obtain all 30 one should cyclically permute the blocks and take 
full inverses of the loops. We give the loops as an edge path as well as their 
decomposition into the basis loops li,l2,h- 

Remark 1.12. It is interesting to note that the explicit isomorphism F3 — >■ 
7ri(r_|-) given on the generators a.j of F3 by a.; H> induces a new length on 
F3 as the minimal length in the letters e^. It is this length that remembers 
the structure of the Gyroid. The computation above is a good example of this. 
The three examples of commutators are minimal elements of [F3,F3] which are 
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loop as an edge path in e^s 



in basis L 



N, 



(hhhr'ihhh) 

[hM] 



eff 



edge missed 



e2(-e4)e5e6(-e2) e3e4(-e6)(-e5)(-e3) 
ei(-e5)e4(-e6)(-ei) e3e5e6(-e4)(-e3) 
ei(-e5)e4(-e6)(-ei) e2(-e4)e5e6(-e2) 
ei(-e5)(-e3) e2{-eQ){-ei) e2,e^ef=,{-e2 
eie6(-e2) e2,ei{-eQ){-ei) e2(-e4)(-e3 
ei(-e5)e4(-e2) e3e5(-ei) e2(-e4)(-e3 

Table 1 . A generating set the shortest loops and their effective 

normal vectors 



(1,-1,0) 
(0,1,-1) 
(1,0,-1) 
(1,1,0) 
(0,-1,-1) 
(-1,0,-1) 



ei missmg 

62 missing 

63 missing 

64 missing 

65 missing 

66 missing 



not the identity in the length in Oj or /j. They are of length 4. In the 6j their 
length is 10. The other examples are of length 6 in the Oj, so not minimal in 
this metric, but they are minimal in the 6j being of length 10. 

1.5.3. Effective normal vector. One distinguishing feature of the different 
loops is their spatial orientation. This for instance has an effect on the flux 
through a surface bounded by such a loop. 

We assume a constant magnetic field. In this case if 5 is a surface bounding a 
loop L = , . . . , ej„ , let fj := '^'l^i 6^^ . Then by Stokes, we can just integrate 
over the surface given by the union of the triangles defined by (/j, /j+i), j = 
1, . . . , n — 1. li Nj = fj X /j+i = fj X Cj+i and N^ff := X^j=i -^j then we have 
for the magnetic flux through S 

ff BdS = Y,lB-N, = ^B.N,ff (1.12) 
JJs J 

The values for N(,ff are listed for the basic loops. Notice that inversion of a 
loop inverts the normal vector, while the cyclic permutation of the components 
leaves the effective normal vector invariant. 

1.6. The quantum graph. In order to promote the skeletal graph to a quan- 
tum graph, we will flx a Hilbert space and a Hamiltonian. Here we follow the 
terminology that a quantum graph is graph with an associated Hilbert space 
and a Hamiltonian on it. In this paper we will use the generalized Harper 
Hamiltonian |2l[20]. 

The original Harper Hamiltonian [12] is obtained for a cubic lattice by using 
the tight-binding approximation and Peierls substitution for the quasi momen- 
tum [16]. 

In the next section, we will give the general theoretical setup for this using 
graphs, groups and representations. For the Gyroid the notions such as graphs 
and groups have been introduced above, so that the reader may substitute these 
in the general definitions below. 

2. Graph Representations and Matrix Harper Operators 

One idea in studying the non-commutative aspects of a given system is to 
give a K-theoretic gap labeling for the Hamiltonians. For this one considers 
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an algebra generated by the Hamiltonian and the symmetries. If everything 
is commutative then this algebra is basically the C* algebra of functions on 
the Brillouin zone/torus. Wc will make these ideas precise using physical and 

mathematical lattices as defined below. 

2.1. Graph language. For us a graph or an abstract graph T will be a collec- 
tion of vertices ^(r) and a collection of edges -E'(r) which run between vertices 
up to bijections preserving the incidences. Each edge can have two orientations. 
An edge together with an orientation is called an oriented edge. Each oriented 
edge e has a starting vertex s(e) and a target vertex t{e). A graph V is called 
finite if both V{V) and £'(r) are finite sets. 

A graph naturally becomes a topological space if the edges are replaced by 
intervals. This space is called the realization of the graph. In more technical 
terms the data above gives a one-dimensional CW complex and we take the 
realization of this complex. When we talk about topological properties of a 
graph, like its fundamental group, we always mean the topological properties 
of its realization. 

A graph is connected if its realization is connected. This means that one 
can travel to all vertices from any given vertex along the edges. A tree is a 
connected graph whose realization is contractible (i.e. the graph has no loops). 
A choice of root of a tree is simply a choice of a vertex and a rooted tree is a 
tree together with a choice of a root. Given a graph F a subgraph r is called a 
spanning tree if it is a tree and the vertices of r arc all of the vertices of F. To 
have a spanning tree F needs to be connected. In this case there are usually 
several choices of spanning trees. A rooted spanning tree is a spanning tree 
together with the additional choice of a root. 

Proposition 2.1. Let T be a finite graph and t be a rooted spanning tree. Let 
vq be the root then 7ri(r) := 7ri(r,t;o) = F„ where F„ is the free group in n 

variables and n = |i^'(r)| — |iiJ(r)|. 

Proof. Consider the graph F/r obtained by contracting all edges of the subgraph 
r. This is the graph which only has one vertex vq and all the edges are loops. 
The two graphs F and F/r arc liomotopy equivalent and hence have the same 
fundamental group. The graph F/r embeds into the plane with n punctures 
with each loop going around one puncture. To obtain a compatible embedding 
of F one blows up the only vertex of F/r into the tree r and both graphs are 
homotopy equivalent to the punctured plane. F is thus homotopy equivalent 
to the wedge product of n circles . It is well known that the first homotopy 
group of this space is the free group in n generators. □ 

An embedded graph is a graph F together with an embedding of it realization 
into an M". Some properties we discuss depend on such an embedding or are 
derived from it. 

Example 2.2. We have so far considered f ^''2'**"' and f + as abstract graphs and 
we have considered F_|_ both as an abstract as well as an embedded graph. The 
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properties of having a certain number of loops at a given point are properties of 
the abstract graph, while the proof made use of the embedding. The effective 
normal vectors are properties of the embedding of the skeletal graph r+ into 

2.2. Graph Harper Operator. 

Definition 2.3. A C* representation pr of a graph F is given by the following 
data. 

• A collection of separable Hilbert spaces one for each v eV{r). 

• A collection of isometries Ug : ^s{e) ~^ -^(e) e^ch oriented edge e 
such that UgUg = 1 whenever e and e' are the two orientations of the 
same edge. 

Remark 2.4. This construction can also be stated in more categorical terms. 
It is a certain quiver representation. A finite graph T also generates a groupoid 
(that is a category in which all morphisms are invertible) and in this setting 
a representation is a functor from the groupoid to the category of separable 
Hilbert spaces. 

Definition 2.5. Given a representation pr of a finite graph V we set := 
®ve.V{f) ^® define the graph Harper operator H on Jiff as 

H:= J2 (2.1) 

oriented edges e 

where we have (ab)used the notation Ug to denote the partial isometry on 
induced by the operator of the same name. 

2.3. Matrix actions. 

2.3.1. General Setup. Fix some finite index set /, a fixed index o e I and 

and order on I, such that o is the smallest clement. Fix isomorphic Hilbert 
spaces i £ I and let (6.; : — > be fixed choice of isomorphisms. We 
allow a choice of 4>o- It might be that (po = id but this is not necessary. 
Set^ = 0.gjJ^. 

In this situation, there is an action of M^j^{End{,Jifo)) on given as follows. 
Let A G M^j^{End{Jifo)) be an endomorphism valued matrix. Then the action 
of A on is given by 

= 4 ^ (2.2) 

Vice-versa given an endomorphism H G End{Jif) there is a corresponding 
matrix M{H) G M|j|(J^) given by 

^ffi|/| ^^^H^®^' ^e|/| (2.3) 

That is if $ = 0^ ^i, M{H) = ^-^oHo^. 
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2.3.2. Matrix action for the graph Harper operator. In order to write 



the Harper operator of ^2.2 as a matrix according to ^2.3.1 we need to fix 
several choices. We will assume that the graph T is connected. First there is 
a choice of base vertex Vq and a choice of order on all vertices in which the 
base vertex is the smallest. Then for each vertex v we choose a fixed path 
Py — that is a sequence of ordered edges ei, . . . , — from Vq to v. We then 
set (j)y := (pgi^ o • • • o : ,J^y^ — )■ J^y. If v = Vq we also allow the choice 
(j)y^ := id : — >• J^vo- This corresponds to the empty path. 

Given such a choice, we obtain the action Pp^^,,,^p„ as above. We call the 
resulting matrix the graph Harper operator in matrix form. 

Given a rooted spanning tree r of F choose vq to be the root and then there is 
a unique shortest path on r from vq to any vertex Vi of r which we can view as 
a path on F. Thus we get the data needed. Again a convenient choice of paths 
is given by a spanning tree. We sometimes write -ffr,T for the corresponding 
graph Harper operator in matrix form. 

Remark 2.6. Any two matrix Harper operators are conjugate. This follows 
simply by making a base change. 

3. C*-GEOMETRY OF HARPER HAMILTONIANS ON LATTICES 

For lattices there are different definitions and connotations in the mathe- 
matics and physics literature. We will use the adjectives "mathematical" and 
"physical" to distinguish the two. As a reference in the physics literature for 
the definition of lattices we use [1] . 

3.1. Translation action. The additive Abelian group acts on itself by 
translation. We denote T^(f) = v + w. For a subset 5 G and a subgroup 
L C we denote by L{S) the set of all translates of points in S under the 
action of all elements of L. 

3.2. Lattices: mathematical, physical and Bravais. 

Definition 3.1. A mathematical lattice with a basis of rank m in M" is an 
injective group homomorphism pi : Z*" ^ M". If m = n then such a lattice is 
called a Bravais lattice with a basis. 

If /O is a mathematical lattice with a basis, we let L := pi{'L'^) be its image 
then this is a discrete subgroup of which is isomorphic to . We will refer 
to L mathematical/Bravais lattice. 

A Bravais subgroup or sublattice is a subgroup L' C L which is a Bravais 
lattice. 

Any Bravais lattice L acts by translations on M" and a primitive cell is a 
fundamental domain for this action. In general a cell is a fundamental domain 
for a Bravais subgroup L' of L. 

Example 3.2. The integer points Z" C are a Bravais lattice which we call 
the crystallographic lattice. A crystallographic cell is a primitive cell for this 
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lattice. The unit cube is a primitive cell. The cube of length 2 is a cell that is 
not primitive. 

As a definition for a lattice in the physical sense, we will take the following 
convention: 

Definition 3.3. A subset A C is a lattice if there is a finite set V C 
and a Bravais lattice L such that A = L{V). 

The tuple (L, V) is called a crystal structure or a lattice with a basi^ We 
will not use the latter terminology as it is mathematically confusing. 

Notice that given A neither L nor V are uniquely defined. We can and will 
mostly choose L to be maximal and V to be minimal or primitive. In this case 
we call (L, V) a primitive crystal structure. Sometimes L is called the space 
group. 

Still then there is of course a choice in fundamental domain and hence a 
choice in V, but for any two such choices there is a unique bijection. In fact 
if we let be a set of orbits of the L action; then it is unique and any V is 
just a choice of representatives. It is also easy to say when L is minimal. L is 
minimal if and only if the number of orbits \ V\ is minimal. 

3.3. Graphs and Lattices. To develop a general theory, we will deal with two 
cases. The usual case is that given a lattice we obtain a graph by adding edges 
to the nearest neighbors. But is is convenient to also allow that we already have 
an embedded graph T\ whose vertices are the set A. In this case, we will let T\ 
stand for this chosen graph. 

3.3.1. Canonical Graph of a Lattice. Given a lattice A without a pre- 
chosen graph define the graph T\ of A to be the graph whose vertices are the 
elements of A and whose edges are the line segments between nearest neighbors. 
This graph is naturally embedded in and we will sometimes make use of this. 

3.3.2. Quotient graphs. Given a choice of symmetry group L for A, we 
also define the quotient graph T/<^(L) as T\/L C M"/L. This is the abstract 
graph whose vertices are given by V and whose edges are given by the set of 
orbits of edges of Fa- 

If L is a maximal group we will just write Fa. 

If the unit cell is a cell for L corresponding to a subgroup L', we also define 
the crystallographic quotient graph of A to be the graph Tj^{L'). 

3.4. Group of a lattice. A lattice A also generates a discrete subgroup of 

as follows: using the natural embedding Fa C we can think of any directed 
edge e of Fa of Fa simply as a vector in M". Then there is an associated 
translation operator Tg. Notice that as vectors the translates of all the oriented 
edges correspond to one another. This means that the set of vectors e is given 
by one vector in each orbit of oriented edges of the embedded graph Fa. This 
set is in bijection with the oriented edges of Fa and we will not distinguish 



'Usually one would include labels such as atom labels to V . 
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between them in the notation. That is we write e for both the oriented edge 
of the abstract graph Fa and the vector in M" it is in bijection with. We wih 
enumerate the set of vectors {e} as el. These vectors then give a symmetric 
group of generators of an Abehan subgroup of M" which we call the lattice group 
T(A) of A. It is the group generated by the ej This group are all the elements 
of the form 

T{A) ■.= {^a,ei\a,eZ} (3.1) 

i 

Proposition 3.4. r(A) is a mathematical lattice. 

Proof. To generate the group we only need one choice of orientation per edge. 
Let ej,j G J be such a choice. This choice defines a map Z"^ — t- M". The image 
of this group is a torsion free Abelian group and hence by the structure theorem 
for Abelian groups a free Abelian group of rank < | J|. □ 

T(A) also naturally acts by translations on M". 

3.5. Fundamental group. The graph T\ is a covering space for T\. Since 



we know the fundamental group of Fa by Proposition (2.1) this information 
can be used to calculate the fundamental group of Fa. In order to do the 
computation, we will have to fix some notation. We fix a rooted spanning 
tree r of Fa and set n := [^^(f a)| — |ii^(r)|. By contracting r we obtain a 



surjection vr : Fa — )• A^^^S"^ which by Proposition (2.1) induces an isomorphism 
on fundamental groups. We fix generators of 7ri(A"5'^) and choose lifts li, . . . ,ln 
of them. By definition each loop k is a sequence of directed edges 6,^ , . . . , ej^.^ .j 

for some k{i). We set k := Ci-^ + ■ ■ ■ + Gi^^^^-y 

Proposition 3.5. Let A. he a lattice with graph Fa and finite quotient graph 
Fa. Let t he a rooted spanning tree for T\ and set n := \E{T\)\ — \E{t)\. If 
the vectors li are linearly independent, then 7ri(FA) = [F„,F„]. Furthermore a 
loop on Fa lifts to a loop on T\ if and only if it traverses each edge the same 
numher of times in each direction. 

Proof. We need to compute which loops on Fa lift to loops on Fa and which 
do not. Given a loop I on Fa it can be written as a word in the basis loops 
li, say w = W - if^^'^ where e(j) = ±1. Fixing a pre-image of vq given by the 
rooted spanning tree, this lifts to an edge path on Fa given by the sequences of 
vectors ej^ , . . . , ej^.^ . Notice that l'^^ gives rise to the sequence —ei^^^-^ , ■ ■ ■ , —eii ■ 
These vectors form a closed loop if and only if they all sum up to 0. Now each 
if^ contributes ib/j to the sum. This means that w lifts to a loop if and only 
if = 0. Since by assumption the li are linearly independent this 

happens if and only if for fixed i the sum of the exponents of the occurring 
li is zero, and this happens precisely if the image of w in the Abelianization 
Ah{¥n) '■= F„/[F„,F„] = is 0. This means that the covering group of the 
cover Fa — )• Fa is Ah{¥n) which is a normal subgroup of F„ and hence the 
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Figure 7. The honeycomb lattice Ahom the generators fi for 
L(Ahon) and the generators of T(Ahon) 

covering group of T\ is [F„,F„]. The last statement follows immediately by 
noticing this is true for the li and hence also for their summands Ci- . □ 

3.6. Main examples. 

3.6.1. The honeycomb lattice. The honeycomb lattice (see Figure [T]) is 
an example of a physical lattice. To make things precise, we let A^ex C 
which is the Bravais lattice generated by — ei := (1,0) and es := ^(1,— \/3). 
We set 62 = — ei — es = ^(l, Vs)- A "dual" Bravais lattice A^^^ generated by 
/2 := 62 — ei = ^{—3, Vs) and /s := 63 — ei = ^(3, \/3) acts via translation on 
Ahex- There are precisely three orbits of this action. These are the A, B and 
C lattices, where we fix the A lattice to be the orbit of (1,0), the B lattice to 
be the orbit of (—1,0), and the C lattice to be the orbit of (0,0). The lattice 
Ahc is then defined to be the union of the A and B sublattices. In A/jc there 
are three nearest neighbors as indicated in Figure [7j 

The symmetry group is embedded as the "dual" Bravais lattice above. 
The group T(A/jc) is again isomorphic to Z^, but it is the original triangular 
lattice generated by ei and 63. See Figure [7| This defines the graph Fa^^^. 

3.6.2. The Gyroid lattice graph. The Gyroid lattice is given by the set of 
vertices V+ of F+ of ^ It is also a physical lattice. Here the symmetry group 
is the space group of I4i32 which is the body centered cubic (bcc) lattice. 
The group T(A) is actually the face centered cubic (fee) lattice as shown in 
Proposition |1.8[ 

3.7. Hilbert space of a lattice. Given a lattice or in general a countable 
set A, we define its Hilbert space to be J^a := /^(A) that is the set of square 
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summable complex sequences indexed by A. There is an alternative way to 
think about a sequence {ax : A G A) as a function on A. Given a sequence as 
above the function is given by /(A) := ax- Vice-versa given ip the sequence is 
obtained by setting ax ■= tp{)^)- 

Given an action of a group G on A there is an induced action of G on /^(A) 
given by gip{X) := 'ip{g^^{X)) ior g £ G and ^p G ^^(A). 



A standard basis for / (A) is given by the functions vi{l') = 6i 



If A is a lattice then there is a unitary representation of A on = /^(A), 
given via the translation operators Ti (f // ) = vi^ii . 

Given a lattice A and fixing a translation group L we can decompose the 
Hilbert space Jif\ by breaking it down in terms of the orbits of L. More precisely, 
labeling each orbit by a vertex in f, we obtain the direct sum decomposition 

^A = 0=^^ (3.2) 

where = l'^{L{v)) for any v which represents the class of v. 

3.8. Partial isometries and projections. In general we get a representation 
of L(A) on Notice that there is in general no representation of T(A) on 
We do have a partial action, which is given by partial isometries. If e is 
an edge from to it; in f then it induces an isomorphism Vg : M'^ — >• simply 
by setting V^{V) = — e). This induces a partial isometry on all of ^(A). 

Another way to describe this partial isometry is as follows. Consider 
Since r(A) D A -after shifting A so that G A- we have that =^^(a) ^ 
and moreover T(A) acts on =^^(a) via translation operators Tg. If P is the 
orthogonal projection of =^^(a) to then Vg = PTgP. If furthermore P^ is 
the projection onto C Hx(l) then we can further decompose the action into 
components U^'"" = PyTgP^. Notice that in principle there could be a directed 
edge e' 7^ e with e' = e as vectors in M". This happens for e.g. for the cube in 



Q.3.1 In this case there will be more components or super selection sectors to 



use physics terminology. 

3.9. Projective representations. Up to now, we have insisted that the trans- 
lation operators are a bona fide representation of the translations groups. 

It turns out that to accommodate such physical data as a magnetic field 
one should only expect a projective representation. This is also in accord with 
general quantum theory, where representations are always only expected to be 
projective. 

This is also one of the sources of non-commutativity. The second being 
that one does not expect the symmetries of the system to commute with the 
Hamiltonian on the nose. One actually has a choice either to preserve the 
commutativity of the symmetry operators or to preserve that the symmetry 
operators commute with the Hamiltonian. We recall some standard facts from 
representation theory. 
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3.9.1. Cocycles. A U (1) 2-cocycle for a group G written as c G Z'^{G, U (1)) 
is a map c: G x G ^ U{1) such that 

c{u, v)c{uv, w) = c{u, vw)c{v, w) (3-3) 

which is called a cocyclc condition. 

Notice that if c is a [/(I) 2-cocycle for G and H C G is a subgroup then its 
restriction c\h : H x H ^ U{1) is a ^7(1) 2-cocycle for H. 

Definition 3.6. A morphism p : G ^ U{J^), is called a projective represen- 
tation with cocycle c if p{1g) = idj(f, that is the identity Iq of G maps to the 
identity operator idj^, and p{u)p{v) = c{u,v)p{u + v). 

Example 3.7. A cocycle c is called trivial if c{u, v) = s{u)s{v)s^^{m)) for some 
group morphism s : G U{1). Given a trivial group cocycle one can always 
perturb or scale an existing representation p to a projective representation pc 
with cocycle c by setting Pc{u) '■= s{u)p{u). With the same formula one can 
scale a projective representation with a cocycle c' to one with the cocycle cc' . 

Remark 3.8. We will be dealing with Abelian groups G, viz. for which 
we will use the usual additive notation of and + for the identity and group 
operation. 

In the case of a free Abelian group and its Hilbert space, given any cocycle, 
one can also twist the standard representation to a projective representation 
with that cocycle. 

Lemma 3.9. Let L ~ Z" C M" 6e a lattice and a G Z'^{L,U{1)). Then the 
operators Ui which operate on M'l via 

Ui{vi>) = a{l,l')vi+v (3.4) 
where vi is the standard basis satisfy 

UiUv = ail, l')Ui+v, UiUv = e{l, l')UvUi with e{l, I') = ^||f^ (3.5) 

Proof. Straightforward calculation. □ 

Lemma 3.10. Let B be a bilinear form on M". Then aB{u, v) := et^("'") is 
a 2-cocycle for R" . If A C M" is a lattice and are generators for this lattice 
then the algebra of operators Ui of the a twisted representation is generated by 
the operators Ui := Ue^. 

Proof. First we calculate 

aB{u,v)aB{u + v,w) = exjp(-[B(u,v) + B{u + v,w)]) 

= eic-p{^[B{u,v) + B{u,w)+B{v,w)]) 
aB{u,v + w)aB{v,w) = exp(-[B(u,v + w) + B{v,w)]) 

= eicp{^[B{u,v) + B{u,w)+B{v,w)]) (3.6) 
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secondly if ^ = ^ OiCi then Yl U^' oc Ui by the formulas (3.5). □ 



3.9.2. Noncommutative Tori. A standard example which will be impor- 
tant in the following is given by the projective action of Z" on J^(Z") with 
the cocycle a given by choosing an anti-symmetric bilinear form B = lirQ 
on M" and then restricting the cocycle. This yields a cocycle which we call 



a. 



e 



aB{u,v) = e^'^®("'^). 



In this case the operators Ui := generate the algebra of operators Ui. 
These generators satisfy 



UiUj = e2"*®»^T,r„ G,,- = G(e„ ej) (3.7) 

This follows from (3.5); e{Ui,Uj) = g^'fQij e"*'^®^* = e^'^*®*J by anti-symmetry 
of 0. In general: 

Definition 3.11. For a fixed n x n anti-symmetric matrix Q the C*-algebra 
generated by n unitary operators Ui on a separable Hilbert space satisfying the 



commutation relations (3.7) is called a non-commutative torus and denoted by 



110. 

Note that T" := Tg is the commutative C*-algebra corresponding to the 
torus T" = (S"^)^" under the Gelfand-Naimark theorem. 

Example 3.12. n = 2: In this case the skew symmetric matrix can be written 
as = ^ . And a{l,l') = exp[iTr6{l A I')], where for I = (Zi,/2) and 

I' = {l'i,l'2)- IM' := det and accordingly e(Z, /') = e^^^^'^''. This case is 

written as T^. 

Remark 3.13. As usual, once a basis hi for is fixed there is a bijection 
between anti-symmetric bilinear forms and skew-symmetric matrices given 
by 0jj = Q{bi,bj). If thus choice of basis has been made, we will write oq. 
In our applications, the basis bi will be given to us by a choice of basis for the 
Br avals lattice L. 

3.9.3. Wannier or Magnetic translation operators. In case of magnetic 
field there is a standard cocycle and representation coming from the i?-field. 
This was first used in A magnetic field B in mathematical terms is a 

2-form on M". 

If we assume that the B field is constant then this is nothing but a skew- 
symmetric bilinear form on M"". Thus giving rise to a cocycle as above. 
This cocycle can now be restricted to any lattice in M". Furthermore in this 
situation, we can choose a magnetic potential A. This is a 1-form on M" 
such that dA = B. This form exists, since dB = and as is contractible 
H*(R'^) = 0, i.e. the reduced cohomology vanishes, so that every closed form is 
exact. 

In this case the cocyle is trivial and the twisted action can be rewritten as 
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Ui,il:{l) = e-^ (3.8) 

Remark 3.14. If X^ILi ~ ^ closed cycle of vectors bounding a simply 
connected polygonal region D with vertices vi, . . . , Vn, s.t. Vi+i = Vi + rui then 
Umi ■ ■ ■ Uran = e~^^id where F is the flux of B through D. Furthermore if B has 
a non-singular vector potential A such that curl{A) = B then F = jj^ BdS = 

nii ■ A{vi + mit) dt 
This means that given two elements 1,1' we have that 

UiUv = e{U)UvUu e(Z,0 = e*^«^''^ (3.9) 
where R is the rectangle spanned by the vectors / and 

Example 3.15. In the square lattice if we choose a constant magnetic field in z 
direction B = k or B = 2Tr9dx A dy then we can choose A = ^{aiy dx + a2X dy) 
with 2719 = a2- ai. Setting U := Ue^ and V = Ue^ we obtain UV = e^'^^^VU 
and this is just e^"^ where 4> = Jd BdA is the flux through the domain full square 
A spanned by ei and 62- The resulting algebra is Tg. 

3.10. Harper operator for a lattice A and with graph Fa. In this section 
we construct a Harper operator for a given lattice A and a graph for the lattice 
Fa. Associated to each such lattice there is a natural separable Hilbert space 
J^\. The Harper operator is an operator on which is obtained by giving a 
graph representation of F. To do this we flx a maximal translational symmetry 
group L for A (see ^ 3.2 ) . Using the techniques of ^we obtain a Harper operator 



as we discuss in detail below. This operator together with the representation 
of the translation group L defines a subalgebra of the endomorphism algebra of 
which we call the Bellissard-Harper algebra and denote by =^(Fa, L, a) or ^ 
for short. Here we allow for pro jective representations acting with the cocycle 



a, e.g. via a magnetic field (see ^3.9.3 ). If F/,, L and a = ob with B = 27ri0 is 
fixed, we write for the corresponding Bellissard-Harper algebra. 

Furthermore again using ^ we can get a a matrix representation for the 
Harper operator and the action of L. These matrices are elements of a matrix 
algebra with coefficients in the algebra generated by the operators corresponding 
to L. Fixing Fa, L, a path basis, a basis for L and a magnetic field B = 27r0, 
we arrive at an algebra of matrices, which we call ^q. The details are given 
below. 

3.10.1. Graph representation of F. To put ourselves in the situation of 



^2.3.1 we will decompose Jif^ as 

^A = 0=^€ (3.10) 

where = l'^{L{v)) for any v which represents the class of v. 

We also fix a 2-cocycle ub corresponding to a skew-symmetric bilinear form 
or equivalently a constant i?--field. This gives us the Hilbert space part of the 
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data of a graph representation pp. We define the isomorphisms as Ug. For this 
we use the fact that an oriented edge of T has a natural representation as a 



vector e (see 3.4) 



Definition 3.16. The Harper operator ii^A is defined to be the graph Harper 
operator Hf on M'\ corresponding to the graph representation pr- 

The Behissard-Harper (BH) algebra =^(rA, -Z^, as) is the C*~algebra of op- 
erators on generated by the projective representation of L and the graph 
Harper operator. 

Remark 3.17. The operator also generates random walks and is related 
to a discrete difference operator as follows. Let Fa be a fc-regular graph, which 
means that each vertex has valence k. Then A = A; — acts as the difference 
operator. A(^)(l) = {Ylie-s{e)=i ^(0 ~ ^(^e(0)) where the sum is over "nearest 
neighbors" as defined by Fa. 

3.10.2. A Matrix representation of the Bellissard Harper algebra. 

In order to get a matrix representation, we fix a vertex Vq of Fa and a choice of 
paths Pv from Vq io v. We will call such a choice a choice of path basis. Again 
a convenient way to fix such data is to specify a spanning tree. 

We then get a matrix representation of the Harper operator and the operators 
coming from the projective representation of L. 

Theorem 3.18. Fixing a choice of path basis, and a basis for L the correspond- 
ing faithful matrix representation of ^{T\, L,aB) is a sub-C* -algebra of 
the C* algebra M|y(f)|(Tg). 

Proof. Before passing to the matrix representation all the operators involved 
are shifted translation operators. Those coming from L and those coming from 
L(A). First we have to show that the operators from L still act as operators 
from L when restricted to J^^, but this is clear since these are diagonal in the 



direct sum decomposition (3.2). Thus the operators in question are conjugates, 



Up^UiU*^ oc Ui for any Ui E Tq. Here G is the matrix obtained from & = 
by using the choice of basis of L. Secondly, for /' G T(A), Up Ui/U*, oc Uii' acts 
as a translation operator which preserves the Vo summand. This means that 
the sum of vectors = —py + V — p'^ is actually in L. Hence the assertion 
follows. □ 

Notice that different choices of path basis may lead to different representa- 
tions, but all these representations are isomorphic; moreover they are conju- 
gates of one another. The effect of changing the basis of L is to replace the 
matrix with its basis transform 0', but as C* -algebras Tq = Tq, — only the 
presentation has changed — with the base change acting as an endomorphism. 

Corollary 3.19. If Q is rational then the spectrum of H\ has finitely many 
gaps. Moreover the maximal number is determined by the entries ofQ. 

Proof. Since there is an injection of into M^y(^p^^{Tq) we can restrict the 
tracial states to The image of the tracial states of Tq is known to be 
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S = 'L + '^ij Oijlj C M ^ We fix a faithful tracial state r and then we have 
that for any gap projection Pa: t{Px) G [0, |T/(r)|] n S. We thus see that there 
are only finitely many possible gaps if all the 9ij are rational. □ 

3.11. Geometry of In general, we are given a lattice A and perhaps the 
graph T\. We can then obtain a family of BH-algebras by choosing different 
cocycles a2^§ • We will call an element of this family . Now we have already 
shown that such an algebra has a faithful matrix representation C MkiT^) 
where k depends on F. It is interesting to note that this family of subalge- 
bras has different geometries and K-theories depending on the choice of 0. 
Generically one would expect 

Expectation 3.20. If Q is generic (i.e. all entries are irrational) then = 
M|y(-p^|(T0) which is Morita equivalent to Tq. 

If this expectation is met is of course dependent on the choices. It is true for 
all the cases we will study. The main motivation is that the non-commutative 
torus at generic G is simple, i.e. there are no two-sided ideals. This usually 
allows one to find that all the elementary matrices are in the algebra and hence 
the algebra is the full matrix ring. The details of our particular calculations 
given in §4 also illuminate this expectation. 

An open question is what happens at non-generic values of 0, i.e. if one or 
more of the entries of B are rational. This again heavily depends on the entries 
of H\. In the cases we study below either = -^[^(riC^e) ^-g^-™ it is a 
genuine subalgebra. This is for instance a good new source of such algebras 
and for families in which the i^T-theory may jump. 

The commutative case is also very interesting. Here we can characterize 
the C*~algebra by the space it represents via the Gelfand-Naimark theorem. 
For this we need some terminology. For each character, or (7*-algebra morphism 
X : T" — 7- C there is an induced C* -algebra morphism x '■ -^/c(T"') — t- Mfc(C) for 
any k. We fix A; = |y(r)|. We say -ffA is generic if it has k distinct Eigenvalues 
in or equivalently if there is a character x of T" such that x{H) has k distinct 
Eigenvalues. 

We call a point x of degenerate if x{H) has Eigenvalues with higher 
multiplicities. The action of L gives an inclusion: i : — )■ ^q. Given a 
character of I^q, we also get a character of T" by pull-back along i. We call a 
point X of ^0 degenerate if i*{x){H) has Eigenvalues with higher multiplicity. 

Theorem 3.21. If H\ is generic then = C*{X) where X is a generically 
k-fold cover of the torus T". This cover is ramified over the locus of degen- 
erate points and is moreover a quotient of the trivial k-fold cover. Here the 
identifications are along degenerate points of X . 

Proof. First notice that the trivial A;-fold cover of T" has the C*-algebra 
T"[ei, . . . , efc]/i? where the Cj are self-adjoint and generate a semi-simple alge- 
bra. This means that the relations R are equivalent to the equations CiCj = SijCi 
and X] 6« = 1- Here the ei can be understood as the projectors to each of the 
copies. 
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Secondly we characterize ^^q. It is certainly a quotient of the C* algebra 
T[H] where H new self-adjoint generator which commutes with the previous 
generators. The way to understand its quotient is as follows. By the theorem of 
Caley-Hamilton we know that the characteristic polynomial p oi H annihilates 
H: p{H) = 0. We claim this is the only relation. Indeed if there were any 
other relation r, then we could write r = r' + r" with r' G {p{H)) and r" a 
polynomial in H of degree less than k. This relation would hold after applying 
any character X) i-e. x{f) = 0- Since = we also get that xi'f") = 0- 

But generically there are k distinct Eigenvalues, so that for generic choices of 
X the minimal polynomial of x{H) is the characteristic polynomial. Thus as a 
polynomial in H the degree of r" must be bigger or equal to k which means 
that r" = 0. Hence = T''[H]/p[H]. 

We can now give the C* morphisms corresponding to the geometric contin- 
uous maps 

UiT" (3.11) 




rpn 

Let Aj € T" be the Eigenvalues of as a matrix with coefficients in the 
commutative ring T*^. Then the map (j) is given by -ff i— >• ^ AjCj. The maps vri 
and Ti2 are just the inclusion maps. There are sections of tti given by sending 
Cj I—)- which give rise to sections = (po Sf^. The corresponding C* map 
is given by i— A^. From this the claims follow readily. □ 

4. Results for the Bravais, Honeycomb and Gyroid Lattices 

4.1. The Bravais lattice cases. 

4.1.1. The case. In case A = Z", we see that L = T{A) = and La is 
the graph with one vertex and n loops. 

From the graph T\ we can already read off that 7ri(rA) = [F„,F„] according 



to Proposition 3.5 since the condition is obviously satisfied. Minimal loops are 
of length 4 and there are (2) unoriented loops. 

Fixing a cocycle a by fixing an anti-symmetric matrix @ (recall that this is 
equivalent to fixing a constant -B-field), the corresponding Harper operator is 
just 

n n 

H^- = Y.Ue, + = Y,U, + U* (4.1) 

i=l 1=1 

where are the standard unit basis vectors of M" and Ui are the generators of 
the non-commutative n-torus Tq. 

The algebra generated by the representation of L = Z" is just the non- 
commutative torus Tq and since H G Tq the algebra is also the non- 
commutative torus ^0 = Tq. 
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Figure 8. The graph T\ for the honeycomb lattice, the span- 
ning tree r (root A, sole edge ei) together with a set of oriented 
edges, the quotient graph Ta/t 



4.1.2. Other Bravais Lattices. If A is the set of points of a Bravais lattice, 
then again L = T(A) = A. For the graph T\ we need the information, which of 
the distances between vertices of A are minimal or we need the additional data 
of a graph Fa. This information is also crucial in determining tti. 

Let ej : j G J he the collection of these minimal vectors and fix an orientation 
ej for each of them. As A = L is a subgroup G A and the minimal vectors are 
given by the A G A with minimal length. Again fix a cocycle a by fixing the 
anti-symmetric matrix 0. Then 

HA = Y,Ue,+Ul^ (4.2) 

The algebra generated by L is always Tq and since H\ E Tq we again obtain 
^0 = Tg. 

Example 4.1. The triangular lattice. This is the lattice spanned by the vectors 
ei = (1,0) and 62 = (^j ■^) in M^. In this case the graph Fa has one vertex 
and three loops with the six oriented edges corresponding to itei, ±62, ±(ei — 



62). Hence the condition of Proposition 3.5 is not met. One can compute the 



fundamental group by elementary methods 
by 



1 

The choice of O is given hy Q = 6 [ ^ ^ ) . The Harper operator is given 



= Ui + Ul + U2 + U^ + + (4.3) 

where U1U2U3 = e^' ^^id. Since H G the BH algebra is = in the 
notation of Example |3.12[ 

4.2. The Honeycomb Lattice. 

4.2.1. Classical geometry. The honeycomb lattice A = Ahon is the two- 



dimensional lattice specified in ^3.6.1 Its quotient graph Fa is the graph with 



two vertices and three edges depicted in Figure pi 
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We set ei = (—1,0), 62 = ^{1,^/5) and 63 = 5(1,— \/3)- We choose the 
rooted spanning tree r (root A, sole edge ei as indicated in Figure [s]). The 
loops a and /? lift as I2 = e2(— ei) and Z3 = e3(— ei), we see that I2 = ^{S, V^) 
and ^3 = 2(3, \/3). Thus the condition of Proposition 



3.5 



is met and we obtain 



Proposition 4.2. 

7ri(rA) = [F2,F2] (4.4) 
At each point of the lattice there are six directed or three undirected minimal 
loops of length six. 

Proof. It is clear by inspection that there are no loops of length two and four 
and the loops through an oriented edge twice will give length bigger than six. 
To get a minimal loop, we need to pass through each edge exactly once in each 
direction, if possible. This can be done by fixing the first edge, and then choos- 
ing among the two left over edges the second returning edge. Now everything 
is fixed. One has to leave the vertex by the last edge and one concludes a cycle 
of 6 edges by choosing the only possible non-traversed edge oriented without 
using an oriented edge twice. This of course gives the known three unoriented 
loops at each vertex. □ 

Remark 4.3. In the generators a, /3, these loops are given by the following 
elements: [h, k]^^ [^2 ^ and [h, 1^^]^^ 

It is interesting to note that the map Fa — )• Fa/t induces a new length 
function on the free group F2 that is not symmetric in symmetric generators. 
E.g. the commutator [^2^"'^, /3^"'^] has length 8 in this metric whereas all other 
commutators listed above have length 6 in the generators Cj. 

Remark 4.4. There is another observation. Given a loop on Fa it decomposes 
into blocks of loops on Fa. Here we could take the loop li = ei(— 62) e3(— ei) e2(— 63) 
which decomposes into 3 blocks. Now cyclically permuting these blocks, we also 
get a loop on Fa and of course the inverse of any loop is a loop, so that the 
loop li generates all six loops. 

4.2.2. Quantum geometry. Fixing a constant magnetic field B = 27r0 
amounts to choosing an anti-symmetric bilinear form Q on M^. We will use 
the induced cocycle both on L and on T(A). The matrix is the matrix of Q 
with respect to the basis (/2, /s) of L. It is completely determined by its value 
0(/2, /s)- The cocycle induced on r(A) will also play a role. It is fixed by the 
value 0(ei, 62). We fix the notation that 

^:=e(/2,/3), q:=e^^'' and = e(-ei, 62), X ■= e'^t (4.5) 
The values of x and 9 are not independent: 

9 = -3(f) and q = (4.6) 

This follows from the elementary calculation 9 = 0(/2,/3) = Q{e2 — 61,63 — 
ei) = B(e2, 63) + 6(62, -61) + Q{-ei, 63) + 6(-6i, -61) = 6(62, -61) - <p + 
e(-ei, -62) = -3(j). 
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The Hilbert space = /^(Ahon) sphts 

'^hon = '^A ffi '^B (4.7) 

where A and B are the vertices as indicated in Figure [8| The six oriented 
edges zbej, i = 1, ... 3 of Fa give rise to three partial isometry operators U±i. 



^' ' c/p, oj' ^-*-lo 



where and U-e^ = U^^ = U*^ are the isomorphisms as in ^3.8 
The Harper Hamiltonian now reads: 



3 



Ue, + + Ues 

i— 1 

In order to put it into a matrix form with coefficients in Tg, we again choose 
the spanning tree r as indicated in Figure |8j 
The Hamiltonian now reads 

HA,r={,^',, ' + '': + neM,(T^) (4.9) 



l+U+V 
where we have used the operators U := xUf2 and V = xUj^- We have that 

UV = qVU or UVU* = qV (4.10) 
The symmetry algebra generated by the translations U is isomorphic to 
on J^. It acts via the representation defined by p(U ) = diag{U , ) 
diag{xU,xU) and p{Uf^) = diag{Uf.^,x^Uf^) = diag{xV,xV). 

Proposition 4.5. Ifq 7^ ±1 then SSq = M2(T|) and hence is Morita equivalent 
to Tl- 

Proof. The strategy of proof is to show that the elementary matrix E12 S SSq- 
In case this happens, we get that all elementary matrices are in since 
E21 = El^, Eu = E12E21 and £^22 = ^21^12 and then SSq = M2(Tl). 

The method to obtain £^12 is by direct calculation using the commutation 
relations (4.10). 

The first step is to set X = p{x^U)p{x^V*)Hp{V)p{U*) and then set Xi = 
[H — qX] . Note we are using the assumption 7^ 1 . 

^ ^ i^-[(i - x') + (1 -t)u* + (1 - x')V' 

In step two, we set X2 = H — Xip{l + ^ + ^)~(q q) jwhere explicitly 
** = BU* + CV* +DU + EV + FU*V + GV*U 
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with 



1_^6 ' l_^6' l_^6 ' l_^6 ' l_^6 

notice that = 9 7^ 1 by assumption. And the coefficients do not vanish if 
7^ 1 and x^ 7^ 1- But x^ = 1 impHes = 1 and x^ = 1 inipUes x^^ = 1 
and so x^ = ^ = So if g 7^ ±1 then all the coefficients are non-zero. 

We can obtain £'12 in several steps by setting X3 = X2 — p{U)X2p{U*), = 
X3 - p{V)Xsp{V*), X5 = p{U)X4piU*) - qXi and ffiially obtain £12 = X5 - 
p{ U*V)X^, where for the last step we need the assumption q 7^ —1. 

The situation for (7 = — 1 is more complicated. Notice that in this case IS 

2 

not simple. For instance there is a * homomorphism of : T 1 — t- Cliff {—Id2)®C 

2 

where Cliff {—Id2) (X" C is the Clifford algebra over C of the standard quadratic 
form given by the negative of the 2x2 identity matrix Id2- If i and j are the 
usual basis vectors then </>([/) = i, (t>{y) = j- 
We let / := kercj) =< 1 + C/^, 1 + V"^ >. 

There is also an algebra involution A on Ti given hy U = —V, V = —V. 



Proposition 4.6. If q = -I and x 7^ 1 then = M2{T{). If q = -I and 

2 

X^ = 1 then SSq is the subalgebra of M2{T\ ) given by matrices of the form 



2 

a b 



b a 



+ J a,beTi,J e M2{^) (4.11) 



Proof. In the case q = —1 we can at first proceed as in the proof of Proposition 
|4.5[ As Xi we obtain 



Here the two cases split dependent on whether x^ = 1 or x"^ 7^ 1- 

We will deal with the case x"^ 7^ 1 first. We set X2 = Xi - p{U)Xip{x^U*), 

X^ = X2+ p{U)X2p{x^U*). Then prr^^^y^jqi^^^s = -£^21 and hence we get 

£'12 = G and hence SSq is the full matrix algebra. 

In case q = —1 and x^ = 1 then since q = x^ = —1, we get x^ = —1 and 

Xi = I = E12 + E21 G ^0. Furthermore, we get that 

'0 l + U* + V* - U -V]^ 
1 



X2. 



Setting yg := \p{U)[X2+ p{U)X2p{-U*)] and = \p{V)[X2- p{U)X2p{-U*) 
we get that 

/O l + C/2\ ^ ~ /O 1 + F2\ 
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Let be the algebra above given in (4.11). It is easy to check that SS' 
is indeed a C* algebra since M^i^^) is a two sided ideal of M2{T1) and A 

2 

commutes with *. Now U* + U and V* +V are both in ^ so that 

H=( l + U* + V*\ ( 0\ , 

\1-U* -V* J^\u + u* + v + v* oj 

Furthermore certainly the operators of L are in so that the inclusion C 
^' holds. 

On the other hand, since L C all the matrices diag{a, a) are in and 
by the above / G so that all matrices of the form ( ^ ) are in ^q. 



Furthermore since the matrices (1 + U'^)E2i and analogously (1 + V'^)E2i are 
in taking products with / = £'12 + £"21; we obtain that all the matrices of 
M2(^) C ^e- Hence that SS' C SSq as claimed. □ 

We now analyze the case of g = 1. Let (x) be the principal (algebra) ideal 
generated by x in Tg. Set Ji = {I + U + V) and J2 = (1 + [/* + F*) and 

</l2 = JlJ2- 

Proposition 4.7. If q = I and x = il ^/^en = C*(X) where X is 
the double cover of the torus x ramified at the points (e^'^^s , e^'^*3 ) and 

If q = 1 and x 7^ ±1 then SSq is equal to the matrix C* -subalgebra o/M2(Tq) 
given by matrices in the C* -subalgebra 



piT^)+[f^ '/J (4.12) 



Proof. The first statement follows from Theorem|3.21 for x = 1- The Eigenval- 



ues of H in Tq are ±^y {1 + U + V){1 + U* + V*), which exist by continuous 
operator calculus as the operator under the square root is self-adjoint operator 
of the form AA* and hence has non-negative real spectrum. The points of the 
space X equivalent to ^0 where these two Eigenvalues coincide are given by 
precisely the points above. 

For X = "Ij we remark that = ^0, by using the involution U 1— >■ 
-U, V ^ -V. 

In case that q = 1, but x =^1, we actually know that x is a sixth root of 
unity and = Cs = e^'^^s or x^ = Cf and in these cases X'* = X 7^ 1- Let 



be the algebra of (4.12). Since = J2 and both of them are algebra ideals 
it is clear that is a C*-subalgebra. Then it is clear that C . In 
order to prove the reverse inclusion, we notice that j^[H - p{U)Hx^p{U*)] = 
{1 + U + V)E2i G ^0. Hence also (1 + [/* + V*)E2i G as well as (1 + 
U + V){1 + U* + V*)En G and {1 + U + V){1 + U* + V*)E22 G ,S§q. This 
together with the action of L shows, SS' C S^q. □ 



In the notation of equation (4.5): 
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Theorem 4.8. If q ^ ±1 the algebra has the Kq group Z © Z. // x = ±1 



then Ko{£^-'^ - 



Proof. For a ^ 1 this directly follows from Proposition 4.5 below. For the 
commutative case q = l,x = =tl this follows from the fact that the double 
cover of the torus obtained by identifying two pairs of points of two tori has 
the corresponding Kq group. □ 

4.3. The Gyroid. 

4.3.1. The lattice and sublattice Hilbert spaces. The lattice F^ has 
the Hilbert space J^r+ = ^^(r+)- 

The subspace J^r+ decomposes naturally into subspaces Jifi where Jifi = 
J^z3{vi) ^■iid 1j^{vi) denotes the set of all translates under of Vi. 

Thus = ®J=o^- This is the composition corresponding to f^''?'''*"'. 

In order to write down the Hamiltonian effectively, we will put together the 
summands in pairs. 

jTr = JTa © JTb © JTc © J^d (4.13) 



J^A = J^(BJ^, J^B = J^i(BJ^, J^c = J^(BJ^5, J^D = J^(BJ^4 (4-14) 

This corresponds to passing from Y'^y^'-"-'' to r+ and to our general setup of 
graph Hamiltonians. 

4.3.2. Cocycles for M"^. As discussed, choosing a magnetic field corresponds 
to a projective representation of F+ by unitary operators on with a cocycle 
a. We recall this here in a more familiar three-dimensional setting. 

A skew symmetric bilinear form B translates to a more well known expression 
as follows. Let i,j, k be the standard unit vectors on and set = B{j, k), 
By := B{k,i) and B^ = B{i,j) then if i? = {Ba;,By,B^) 

B{m, m') = B ■ (m x m'), ^^(m, m') = exp{i-B ■ (m x m')) (4.15) 
Indeed is a cocycle since 

B ■ (m X m') + B • ((m + m') x m") = B ■ (m x m' + m x m" + m x m") 

= B ■ {m X [m! + m")) + B ■ {m! x m") 

The physics interpretation of this is that B \s a, fixed constant magnetic field 
and then aB{m,m') is the magnetic flux though the triangle spanned by m an 
m' . 
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4.3.3. Cocycles for the Gyroid. The maximal translation group L for 



is the bcc lattice spanned by the vectors fi given in (1.4) ot gi ( 1.5 ). The lattice 
group T'(r+) is the fee lattice spanned by the vectors 64, 65, eg- By restricting a 
cocycle given by a vector on M^, we obtain a cocycle for each of these lattices. 
Since both the fee and the bcc lattices are Bravais lattices, these are precisely 
all cocycles coming from cocycles based on anti-symmetric forms Q = i^B 
stemming from a magnetic field B. 

We will use the basis {gt) to fix the matrix representation B of 0. 

4.3.4. Graph representation and Graph Harper Operator on ■^t+- 
We fix a vector B. This fixes the corresponding cocycle a{v,w) := e*2^'(''^'^). 
Now consider the graph representation as defined in §3.10[ where the order of 
the vertices is A, B, C, D. 

Using the fixed cocycle a above, we obtain the partial isometries correspond- 
ing to elements of T{T^). These are Ui := C/gi according to the list (1.6) as 
discussed in 



In Hilbert space decomposition (4.13) the Graph Harper Operator Hp be- 
comes the 4x4 matrix. 



Hr 



/o 

U2 

-1 _ 



'1 




TI* 
^6 



(4.16) 



recall that Ui = p{ei) and U* = 

4.3.5. Matrix Harper operator. We choose the rooted spanning tree r 
(root A, edges 61,62,63) as indicated in Figur e [s) U sing this we obtain the 
following matrix Harper operator according to §3.10.2 



/O 1 1 1 \ 

1 U*,UiU2 U^U^Us 

1 U^UeUi U^UiUs 

\l U^U^Ui UIUIU2 / 



/o 
1 
1 



1 



A* 
B 



1 

A 


C* 



1\ 

B* 

C 
0/ 



(4.17) 



The coefficients can be expressed in terms of the operators of the magnetic 
translation operators of the bcc lattice. We fix U := Uf^^V := Uj^ and W := 



C//3 for the fi listed in (1.4). Then 
A = aV*W, B 
where a = "(^^.-e^^Ce^-e^.-eO ^ ^ 

a{-f2,h) ' 



= bWU\ C 



= cW*UV 
and c 



(4.18) 

a(-e2,e4)a(-e2+e4,e3) 
a(-/3,/l)a(-/3+/i,/2)- 



4.3.6. Choices and Notation. In order to proceed we fix some convenient 
basis and notation. First we base change from the basis fi to the basis 

5l = -/2 + /3, 52 = -/l + /3, 53 = /l + /2-/3 (4.19) 



THE GEOMETRY OF THE DOUBLE GYROID WIRE NETWORK 35 

for the bcc lattice B. We see that A = a'Ug^,B = h'Ug^^C = c'Ug^ again for 
fixed constants a' = a(— ei, —ee)a{—ei — ee, 62), 6' = «(— 63, — e5)a;(— 63 — 65, 62) 
and c' = — a( — 62, e/i)a{~e2 + 64, 63). The operators A, B, C again generate the 
C* algebra of the non-commutative torus Tq where Q is the matrix of the 
bilinear form in the basis gi. Explicitly: 

012 = ^B ■ {gi X 52), ^13 = ■ {gi x 53), ^23 = • {92 x 53) (4.20) 

ZVr ZTT ZTT 

We will fix the notation 
ai := e^"^- = a\g,,g2),a2 := e^-'^^^ = a\gi,g^),a^ := e^-^^^ = a'^g^^g^) 

(4.21) 

This means for the commutators 

AB = aiBA, AC = a2CA, BC = a^CB (4.22) 

4.3.7. Matrix action of r(r+). Given the choice of the spanning tree r 

also determines the action of T{T-^-). 

It will be convenient to introduce forth roots of the af. 

</,l = ef^^l^ (t)2 = e^'^^\ 03 = et^^23^ $ = </,i<^2<^3 (4.23) 
and the matrices 

Ai = diag{l, 01, 02, 0102), A2 = diag{l, 0i, 0i03, 03), A3 = diag{l, 0203, 02, 03) 

(4.24) 

By definition 0:10203 = 

Then the action p is given by the following matrices 

p{A) = Aidiag{A, A, A, A),p{B) = A2diag{B, B, B, B), p{C) = K^diag{C, C, C, C) 

(4.25) 

We make the following observation that makes the calculations quite a bit 
simpler. 

AiA^ = dia5(l,l,$,$), 1^2^-1 = diag{\,^,^,\), kiKl = diag{l,^,l,^) 

(4.26) 

4.3.8. Calculation of and K^{S8q). . 

Proposition 4.9. If ^ ^ 1 then the BH-algebra SSq is the full matrix algebra 

= M4(T|). 

Proof. The strategy is to again show that enough elementary matrices are in 
This will be done in a case by case study. We will present the first case in 
detail. 
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Assume that ai / <I>^, 02 / ^* and 03 / Then wc have to do 6 steps to 
obtain £'34. These are 



The resulting matrix is 

(1 - ^^a2a3)(l - <^>a3)(l - ^a2){^ai - l>^a2a3)($Q!3 - ^as){^a2 - ^a2)E4s 

and the factor is invertible by assumption. This provides £^43, £^34 = £43, £33 = 
E34E4S and E44 = E43E^. Now to get the other elementary matrices, we first 
assume (pi ^ 1 then set Yi = HE 44- (j)zp{C)E:i4 and Y2 = cp^Yi- p{B*)Yip{B) = 
4>'i{4'i ~ l)-£'24) and Ei4 = Yi — (j)^p{B)E24 which guarantees that all the Eij G 
^e- If = 1 then we use Z'^ = 4)^Yi - p{C)Hp{C*) = ^3(1 - az^E24, so 
since 0:3 7^ $ we obtain E24 and £14 (as above) and hence again the whole 
matrix algebra. The arguments for the cases ai 7^ $,0:2 7^ ^^,ck3 7^ ^ and 
ai 7^ 0:2 7^ ^1 Q;3 7^ are analogous. 

Now by pure logic it follows that is the full matrix in all cases except the 
three cases = aj = = ^ for {i, j, k} = {1, 2, 3}. We treat only one of 
these, as the other two follow symmetrically. So assume ai = 02 = 03 = 
In this case the story is similar, but faster. 

We keep X_i,X2 as above. Then if $ 7^ $ we set X'^ = ^-p{BC*)X2p{CB*) = 
(1 — <I>)(1 — — ^>)£i3 and again obtain the full matrix algebra using HEn 
to obtain the needed elementary matrices. If $ = $ then since $ 7^ 1 we have 
$ = — 1 and hence a2 = = — 1 and so (f)2 7^ 0^ since else 02 = il and 
as = 0t = 1. Thus we get that X^ = ^2X2 - piA)X2p{A*) = 4(<^2 - <^2)^i3 
with the factor being invertible. And again, we are done. □ 

Proposition 4.10. // $ = 1 and not all ai = ai = ±1 then the BH-algebra 

^0 is equal to the full m,atrix algebra = M4(T|,). 

Proof. As before this follows from an explicit representation of the generators 
Eij as expressions in ^q. Since not all ctj = 1 there must be some 7^ 1. 
We will assume that 03 7^ 1. The other calculations arc symmetric. Again one 
has two cases. Either 02 = ^2 or not. In the second case, we set Yi = H — 
p{BC*)Hp{CB*) and Y2_= a2>l = p{AC*)YipiAB*).We set Y3 = a3^i0i^2 - 
p{C)Hp{C*), and Y4 = ^1 021^3 - p{A)Hp{A*) = (1 - a3)(a3 - a3)(a3<^20i - 
0|03)(0i02 -0(2- (l)2(pl)E43. Now as Q!3 7^ 1: Y4j^0 if (1)2 7^ 1- In case 02 = 1, 
we must have 0i 7^ 1 because otherwise 03 would also be 1. In this case, we 
proceed analogously to obtain E24. In both subcases we obtain the full matrix 
algebra following the strategy used in the previous proof. 
Finally, we deal with the case that a2 = 02: 



^1 
X2 
Xs 
X4 
X5 



H - p{AB*)Hp{BA*) 
Xi - p{BC*)Xip{CB*) 

X2-p{AC*)X2p{CA*) 



$01X3 - p{AB*)X3p{BA*) 
$03^4 - p{BC*)X4p{CB*) 
^(X2X^- p{AC*)X^p{CA*) 



THE GEOMETRY OF THE DOUBLE GYROID WIRE NETWORK 



37 



If a2 = 1 then qi = Q3 then it follows that ai / ai and we are done by an 
analogous calculation. Indeed if ai = ai then also 0:3 = 0:3 and all ai are real 
which we excluded. 

If a2 = — 1 then ai = —03 and it follows that ai 7^ cti and we are done by 
an analogous calculation. Indeed if ai = —as then ai = ai again means that 
all three ai are real which was excluded. □ 

Proposition 4.11. //$ = 1 and all of the a-i = 1 then the BH-algebra is 

a proper subalgebra of the full matrix algebra SSq C M4(T^). 

Proof. We will show this using the character % : — > C, defined by xi^) = 
x{B) = x{C) = 1. This character induces an algebra morphism x '■ M4{T^) 
M4(C) and we will show that xi^e) S M4(C) which implies the result. We 
note that 



Since a,- 



V 

all the 



1\ 

1 

1 

0/ 



F - Id, with F 



/I 
1 
1 

V 



1\ 

1 

1 



must be fourth roots of unity and furthermore 
4>i4>24>^ = 1- There are again three cases. 

The first is that all (f)i = l. In this case x(p(^)) = x{p{B)) = x{p{C)) = Id- 
Therefore xi^e) =< x{H) >C M4(C). Now = 4""!^ and so x{HT = 
-^<^+ (EILi so *^at x{^e) is the 2-dimensional subalgebra of M4(C) 

spanned by Id and F. 

The second case is that all 0f = 1 and only one 0i = 1. We will assume that 
01 = 1, (/'2 = ^1; (f's = ~1- The rest of the cases arc symmetric. 

In this case the image of x is the 6-dimensional subalgebra of matrices of the 
form 



(a 
b 
d 

\d 



b 
a 
d 
d 



c 
c 
e 
f 



c\ 

c 

f 



(4.27) 



Let be the sub-algebra above. It is an exercise to check that this is indeed 
a subalgebra. Then as x(p(^)) = x(p(-S)) = xip{C)) = diag{l,l, —1, —1), so 
that xi^e) C On the other hand Pi = EU+E22 = \{Id+diag{l, 1, -1, -1)) G 

and P2 = Ess + E44 = \{Id- diag{l, 1, -1, -1)) e xi^e)- And hence all 
the Pix{H)Pj,i,j = 1,2 G x('^e)- These are the 2x2 block matrices making 
up x{H)- But the span of these matrices is precisely 

The final case is exactly one = 1. We will treat the case 0i = 1,^2 = 
i, (j)s = —i- The rest are again symmetric. In this case xi^o) is the 10- 
dimensional matrix algebra of matrices of the form 



( a b c d\ 

b a c d 

e e g h 

\f f k IJ 



(4.28) 
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Let SS" be the subalgebra above. Again it is an exercise to check that SS" is 
indeed a subalgebra. In the current case: x{p{A)) = x{p{B)) = x{p{C)) = 
diag{l,l,i,-i) and hence it follows that xi-^e) C . Now x(/>(^))^ = 
diag{l, 1, —1, —1) so that the Pi, P2 above are in x{^&)- But furthermore we 
have that ^{x{p{A))P2x{H)P2-iP2x{H)P2 = and ^^{x(p{A))P2X{H)P2 = 
iP2x{H)P2 = -E'43. So that -E33, £"34, -£43, £^44 G x{^b)- The other generators 
of now are given by A, Pix{H)Pi, P2X{H)Pi, EuP2X{H)Pi, PixiH)P2E33 
and Pix(if)P2-E'33 which are all in x{^q)- 

□ 

The trickiest case is the case in which all the CKj are real and not all equal to 
one. 



Proposition 4.12. Let $ 

all the (pi are different 



1 and all the ai real, hut not all ai = 1. Then if 
M4(Te), else C M^{Jq). 



Proof. Since not all ai = 1 there must be exactly one a.; = 1 with the other 
two being equal to —1. We will deal with the case 03 = l,a2 = as = — 1. 
The others are similar. Consider X = H + p{BC*)H p{CB*). X'^ = l2Eii 
and X2 = HEii = E21 + E-^i + E^. We know that 4>i / 4>2 and (pi 7^ (ps and 
X3 := (P1X2 - p{A)X2p{A*) = {(pi - (p2)Esi + {(pi - (p3)E^i. Here we get the two 
cases. If (p2 7^ 03 then we obtain (p2Xs — p{A)Xsp{A*) = {(p2 — (p3){4>i — (p3)E3i 
and we have successively £31, £41, £21 and their transposes in ^q. But this is 
a set of generators. 

In case 02 = 03, we see that Pi = ^^^^^ X2 = £31 + -£41. 
£12, £21, £22 G ^0. Furthermore P2 = Pi £12 = £32 + E42 e 
Id - En - E22 = E33 + Eu e^e,Q2 = QiHQi. 

Cliff {(^ °' 



We shall use the morphism % : Tq — >■ CI 



Hence also 
?0 and Qi = 

C given by 



^0 1^ 

x{A) = x(-B) = ei, x(C) = e2 where ei,i = 1,2 are the generators of the 
Clifford algebra which satisfy ef = \ and eie2 + 6261 = 1. The map x induces 
an algebra morphism x '■ ^^{Tq) ^ CI. 



Claim: The image of x is given matrices of the form 

f a b e e\ 

c d f f 



9 
\9 



h 
h 



(4.29) 



kj 



which is a free rank 10 CI proper submodule of M/i[Cl) and in turn is a 
proper submodule of M4(T|). 

To prove the Claim, we again check both inclusions. Let be the subalgebra 
above. Notice that Id G ^' 



X{H) 



1^ 
1 

1 



1 



ei 
ei 



1 

ei 


62 



1\ 

ei 

62 

0/ 



G 
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XP{A) = XP{B) = eidiag{l, (/>!, 02, (/>2) G ^' and xp{C) = e2diag{l, (/>i, (/)2, 4>2) G 
So that xi^e) C S^'. On the other hand the algebra is generated 
by En,Ei2, E21, E22, Pi, Pf , P2, P2 , Qi,e2XiQ2) = E34 + E43 which are ah in 
xi^e), so that C xi^e)- □ 

Theorem 4.13. If ^ ^ 1 or ^ = 1 and at least one 7^ 1 and all (pi are 

different then = M4(T|) and K{.^e) = Z^. 

If (pi = ^ for all i then K{^q) = K{X) where X is the cover of the three 
torus given by Theorem 3.21\ 



Proof. The only thing that remains to be proved is that H is generic. Indeed 
this can be done by direct computation. We will not give the Eigenvalues here 
as they are quite long expressions. But notice that for x{^) = ~^iX{B) = 
^,x{C) = —1 the matrix x{H) has the four distinct Eigenvalues =b\/5, ±1. □ 

Acknowledgments 

We would like to thank J. Bellissard, M. Dadarlat, M. Kontsevich, M. Mar- 
colli and A. Schwarz for enlightening discussions. Foremost, we wish to thank 
H.W. Hillhouse for sparking the interest in this project and many discussions. 

RK thankfully acknowledges support from NSF DMS-0805881 and the Hum- 
boldt Foundation. He also thanks the Institut des Hautes Etudes Scientifiques, 
the Max-Planck~Institute for Mathematics in Bonn and the Institute for Ad- 
vanced Study for their support and the University of Hamburg for its hospital- 
ity. While at the IAS RK's work was supported by the NSF under agreement 
DMS-0635607. 

BK thanks the Department of Physics at Princeton University for its hos- 
pitality and thankfully acknowledges support from the NSF under the grant 
PHY-0969689. 

Any opinions, findings and conclusions or recommendations expressed in this 
material are those of the authors and do not necessarily reflect the views of the 
National Science Foundation. 

References 

[1] N. W. Ashcroft and N. D. Mermin Solid State Physics, Brooks/Cole 1976 

[2] J. Bellissard, Gap labelling theorems for Schrodinger operators, in: From number theory 

to physics (1992) 538-630 
[3] J. Belhssard, A. van Elst,H. Schulz-Baldes. The noncommutative geometry of the quantum 

Hall effect. Topology and physics. J. Math. Phys. 35 (1994), 10, 5373-5451 
[4] A. L. Carey, K. C. Hannabuss and P. McCann, Quantum Hall Effect on the Hyperbolic 

plane, Commun. Math. Phys. 190 (1998) 629-673 
[5] A. Connes, Noncommutative Geometry. Academic Press, Inc., San Diego, CA, 1994. 
[6] G.A. Elliott. On the K -theory of the C* -algebra generated by a projective representation 

of a torsion-free discrete abelian group. Operator algebras and group representations, Vol. 

I (Neptun, 1980), 157-184, Monogr. Stud. Math., 17, Pitman, Boston, MA, 1984. 
[7] G. A. Elliott and H. Li, Strong Morita equivalence of higher-dimensional noncommutative 

ton. U, Math. Ann. 341 (2008) 825-844 
[8] T. Hahn (editor) International Tables for Grystallography, Vol. A, Springer, 2005 



40 RALPH M. KAUFMANN, SERGEI KHLEBNIKOV, AND BIRGIT KAUFMANN 

[9] K. Grofic-Brauckmann and M. Wohlgemuth, The Gyroid ts embedded and has constant 
mean curvature companions^ Calc. Var. 4 (1996) 499-523 
[10] K. Grosse-Brauckmann, On Gyroid Interfaces, J. Colloid and Interface Science 187 (1997) 
418-428 

[11] D. A. Hajduk ct. al., The gyroid - a new equilibrium morphology in weakly segregated 

diblock copolymers Macromolccules 27 (1994) 4063-4075. 
[12] P.G. Harper. Single band motions of conduction electrons in a uniform magnetic field.. 

Proc. Phys. See. London A68 (1955), 874-878 
[13] S. Khlebnikov and H.W. Hillhouse. Electronic structure of double- Gyroid nanostructured 

semiconductors: Perspectives for carrier multiplication solar cells Phys. Rev. B 80, 11 

(2009) 115316 

[14] C. A. Lambert, L. H. Radzilowski, E. L. Thomas, Triply Periodic Level Surfaces as Models 
for Cubic Tricontinous Block Copolymer Morphologies, Philos. Transactions: Mathemat- 
ical, Physical and Engineering Sciences, Vol. 354, No 1715, Curved Surfaces in Chemical 
Structure (1996) 2009-2023 

[15] M. Marcolli and V. Mathai, Towards the fractional quantum Hall effect: a noncommu- 
tative geometry perspective, in: Noncommutative geometry and number theory: where 
arithmetic meets geometry, Caterina Consani, Matilde Marcolli (Eds.), Vieweg, Wies- 
baden (2006) 235-263 

[16] G. Panati, H. Spolm and S. Tcufcl, Effective dynamics for Block electrons: Peierls sub- 
stitution and beyond. Comm. Math. Phys. 242 (2003) 547-578 

[17] M. Pimsner and D. Voiculescu. Exact sequences for K-groups and Ext-groups of certain 
cross-product C* -algebras. J. Operator Th. 4 (1980) 93-118. 

[18] M. A. Rieffel, C* -algebras associated with irrational rotations, Pac. J. of Math. 93 (1981) 
415-429 

[19] M. A. Rieffel and A. Schwarz, Morita equivalence of multidimensional noncommutative 

tori. Int. J. Math. 10 (1999) 289-299 
[20] T. Sunada, A discrete analogue of periodic magnetic Schrodinger operators, Contemp. 

Math. 173 (1994) 283-299 
[21] A.H. Schoen. Infinite periodic minimal surfaces without selfintersections. NASA TN- 

D5541 (1970). 

[22] M. F. Schulz, et. al.. Epitaxial Relationship for Hexagonal-to-Cubic Phase Transition in 

a Book Copolymer Mixture. Phys. Rev. Lett. 73 (1994) 86-89. 
[23] V.N. Urade, T.C. Wei, M.P. Tate and H.W. Hillhouse. Nano fabrication of double- Gyroid 

thin films. Chem. Mat. 19, 4 (2007) 768-777 
E-mail address: rkauf manOmath . purdue . edu 

Department of Mathematics, Purdue University, West Lafayette, IN 47907 and 
School of Mathematics, Institute for Advanced Study, Princeton NJ, 08540 

E-mail address: skhleb@physics . purdue . edu 

Department of Physics, Purdue University, West Lafayette, IN 47907 
E-mail address: ebkaufma@math.purdue.edu 

Department of Mathematics and Department of Physics, Purdue University, 
West Lafayette, IN 47907 and Department of Physics, Princeton University, 
Princeton NJ, 08544 



